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We report on the formation of surface insta- 
bilities in a layer of thermoreversible ferrogel 
when exposed to a vertical magnetic field. Both 
static and time dependent magnetic fields are 
employed. Under variations of temperature, the 
viscoelastic properties of our soft magnetic mat- 
ter can be tuned. Stress relaxation experiments 
unveil a stretched exponential scaling of the 
shear modulus, with an exponent of /3 = 1/3. 
The resulting magnetic threshold for the forma- 
tion of Rosensweig-cusps is measured for differ- 
ent temperatures, and compared with theoreti- 
cal predictions by Bohlius et. al. in J. Phys.: 
Condens. Matter., 2006, 18, 2671-2684. 



1 Introduction 

Ferrofluids have found widespread applications 
in rotary seals, loud speakers, or medicine [1, 
2]. They are a colloidal dispersion of magnetic 
nano particles in a carrier fluid like water or 
kerosene [3]. Due to their superparamagnetic 
behaviour, they can be controlled by externally 
applied magnetic fields. 

The same is true for ferrogels [4-7], with the 
difference that here the magnetic nano particles 
are embedded in an elastic (polymer) matrix. 
This is handy for many upcoming applications. 
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such as soft actuators, magnetic valves, magne- 
toelastic mobile robots [8], artificial muscles [9] 
or magnetic controlled drug delivery [10-13]. 
Most of the applications are based on the con- 
tinuous deformation of ferrogel bodies in gradi- 
ent fields [14] or in homogeneous fields [15,16]. 
Another promising effect, which so far has not 
been exploited, is the drastic shape transition 
of ferrogels in the succession of the Rosensweig 
instability. If a threshold of the magnetic in- 
duction is surpassed, the fiat reference state of 
the ferrogel is predicted to develop surface pro- 
tuberances, pointing in the direction of the ver- 
tically oriented field. The critical induction B,, 
is derived by means of a linear stability analysis 
in Ref. [17]. It depends on the surface tension 
cr, the mass density p, the gravity Qq, and on the 
elastic shear modulus G of the gel according to 
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where /Xq is the vacuum permeability and /i^ is 
the relative permeability of the ferrogel. In com- 
parison to the Rosensweig instability in ferroflu- 
ids {G = 0) [18] the threshold Be is enhanced by 
the elasticity due to an increase of the surface 
stiffness. In contrast, the critical wavelength is 
independent of the elastic modulus of the gel 
and is described by the capillary wavelength 
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familiar from ferrofluids [18]. Like in the case of 
ferrofluids, neither the critical field nor the crit- 
ical wavelength depend on the viscosity. Note 
that, in contrast to ferrogels, for the instability 
of a ferrofiuid layer covered with a thin elastic 
film, it was predicted, that kc decays drastically 
with increasing elasticity [19]. 
More recently, the authors of Ref. [17] derived 
the final pattern which forms in ferrogels [20]. 
These results are limited to patterns of small 
amplitude and were obtained via a minimisation 
of the energy density, similarly to the method 
used for ferrofluids [21,22]. For ferrofluids, the 
energy density comprises hydrostatic, magnetic 
and surface terms. To account for the elastic 
deformations in ferrogels, an energy contribu- 
tion as given in [23] is added. The Rosensweig 
instability in standard ferrofluids is associated 
with a transcritical bifurcation and exhibits a 
hysteresis [21,24]. In ferrogels, this hysteresis is 
expected to shrink with increasing shear modu- 
lus [20]. 

The Rosensweig instability in ferrofluids has 
been studied in many experiments; see e.g. 
Refs. [18,25-31]. However, its counterpart in 
ferrogels is still awaiting measurements. The 
reason for this void is that up to now mostly 
covalent cross linked polymer gels [5,7,32] have 
been synthesised. This process results usually 
in rather "hard" gels. Due to their high elas- 
ticity and the saturation of the magnetisation, 
one can not excite surface instabilities in these 
gels even for very high magnetic field strength. 
Only recently soft ferrogels have been created 
[6] which take advantage of thermoreversible, 
i.e. physically crosslinked gelators. Contrary 
to chemically crosslinked, irreversible ferrogels, 
their elasticity can be controlled by a thermore- 
versible sol-gel transition. In the following we 
investigate the Rosensweig instability in such a 
thermoreversible ferrogel. 

2 Material and Methods 
2.1 Synthesis 

We prepared a thermoreversible ferrogel by 
swelling 5wt.% of a gelator in an oil-based fer- 



rofiuid containing 30wt.% of magnetite parti- 
cles. The carrier liquid for the preparation of 
the ferrofiuid was paraffin oil (Finavestan A50B 
from Total Deutschland GmbH) with a kine- 
matic viscosity u = 13.6mm^/s at 20 °C and 
a molar mass of 280g/mol (manufacturer in- 
formation). The magnetite particles were pre- 
pared by co-precipitation and stabilised with 
oleic acid [6]. Transmission electron micro- 
graphs show, that the diameter of the particles 
is 8 ± 1 nm [33] . As a gelator we have utilised 
Kraton G 1726 from Kraton Polymers, Belgium, 
which is a mixture of 30wt.% triblock copoly- 
mer and 70wt.% of a diblock copolymer. The 
triblock copolymer is poly(styrene-b-(ethylene- 
co-butylene)-b-styrene) (SEES) with a molar 
mass of Mw = 77 700 and a polydispersity 
index PDI = 1.01 [33]. The diblock copoly- 
mer (SEB) is exactly one half of the triblock. 
The size distribution has been established by 
means of size exclusion chromatography. For 
both the di- and triblock copolymer the styrene 
content amounts to 30wt.% (manufacturer in- 
formation). In comparison with the pure tri- 
block copolymers used in earher studies [33] , the 
diblock acts here plasticiser and lowers the 
softening temperature to 25 °C, according to the 
falling ball method [6]. Both the ferrofiuid and 
the ferrogel reveal no structure in optical mi- 
crographs, i.e. they are perfectly homogeneous 
down to a sub-micrometer scale. The ferrogel 
sample remains stable for 1^/2 year without any 
separation of the fluid phase from the gelator, 
i.e. without any syneresis. Likewise, we did not 
observe changes of the magnetic properties and 
the microstructure within this time. 

2.2 Material properties 

For the characterisation of the elastic proper- 
ties we utilise a rheometer (MCR301, Anton 
Paar) in cone-and-plate geometry. The cone has 
a diameter of 50 mm and a base angle of 1°. 
The rheometer is equipped with a Peltier ther- 
mostated temperature device (C-PTD 200/E). 
Figure [1] displays the result of oscillatory mea- 
surements of the shear modulus at 1 Hz for a de- 
formation of 7 = 0.01 versus the temperature. 
G' and G", which denote the real and imag- 
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Figure 1: Storage modulus G' (red circles) and 
loss modulus G" (blue squares) versus the tem- 
perature for the thermoreversible ferrogel. 



inary part of the shear modulus, respectively, 
have a crossover around 31 °C. Above ~ 45 °C, 
the sample becomes liquid. 

To measure the magnetisation curve of the fer- 
rogel, we produce a sphere by casting 1.0 ml of 
the liquefied gel into a hollow aluminium mold 
at 55 °C. For details see Ref. [16]. The magneti- 
sation of the spherical sample was then mea- 
sured by means of a fluxmetric magnetometer 
(Lake shore, model 480) at ^ = 20 °C. Figure [2] 
shows this data. The sample is superparamag- 
netic with an initial susceptibility of Xo = 0.82. 
This data has been fitted with the model of 
Ivanov [34], assuming a gamma distribution of 
the particle diameters. We use this model as 
well to extrapolate M{H) for all sample tem- 
peratures, by evaluating the fit with a different 
6, while all other parameters are held constant. 
As an example, the dashed line gives the extra- 
polation for the maximal applied temperature 
of ^ = 38 °C. 

The density p of the ferrogel is measured by im- 
mersing the sphere in water with added salt. 
The salt establishes a vertical density gradient 
which keeps the ball floating in the middle of 
the container at half depth. With the help of a 
syringe, water has been sampled directly above 
and below the ball. The density of both samples 
is then measured by means of a vibrating tube 
density meter (DMA 4100, Anton Paar Co.). 
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Figure 2: Magnetisation as a function of the 
applied magnetic induction for the thermore- 
versible ferrogel. The solid line represents a 
least squares approximation with the model put 
forward by Ivanov [34]. The best fit parame- 
ters for the gamma distribution are: exponent 
a = 8.1 ± 1.3, mean diameter xq ■ {a + 1) = 
8.1 ± 0.4 nm, volume fraction = 3.9 ± 0.09 %. 
The core magnetisation has been held constant 
at Mp = 446 kA/m. The dashed line indicates 
the prediction for 6^ = 38 °C with the above pa- 
rameters. 
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From the mean of both values follows the den- 
sity of the ball to be p = 1085 ± 1 kg/m^. 
More cumbersome is an estimate for the surface 
tension a. As a rough estimate, we measure the 
surface tension of the paraffin based ferrofluid 
the gel was created from with a ring-tensiometer 
(LAUDA TE 2). We get ctff = 28.7 mN/m 
for both the ferrofluid and also the underly- 
ing paraffin oil that was used as a carrier liq- 
uid for the ferrofluid. For obvious reasons, the 
tensiometer can not be directly applied to the 
ferrogel. Because the gelator is not surface ac- 
tive, it shall not have a signiflcant influence on 
o. We therefore use ctff in the subsequent cal- 
culations as the surface tension of the gel. 

2.3 Setup for measurements of 
the normal field instability 

Our experimental setup is shown in Fig. [3](a). 
In the centre is a thermostated aluminium block 
of size (x,y,z) = (60,60,50) mm^. A cylindrical 
bore of 50 mm and a depth of 25 mm serves 
as a vessel for the ferrogel. A bore with the 
same diameter, but a depth of 23 mm pene- 
trates the block from the lower side, in this 
way leaving a bottom plate with a thickness of 
2 mm (cf. Fig. [3]a). Utilising thermal grease, 
the vessel is thermally connected to two Peltier 
elements, as shown in Fig. [3](b). They are 
equipped with heat exchangers (lA cooling Co., 
type 1A-SL2), which are circulated by water 
from a closed cooling system (LAUDA RK20 
KP). A thermo-resistor PtlOO serves to moni- 
tor the temperature of the vessel. The Peltier 
elements are powered by a DC-current source 
(EUROTEST Co., type LAB/SL 230/AI/LT) 
which is controlled via IEEE from the computer. 
By a proportional-integral method, the com- 
puter regulates the temperature of the vessel 
with a precision of 5 mK in the range of —35 °C 
to 110 °C. The vessel is covered from above by 
means of an aluminium plate with a thickness 
of 0.5 mm. This lid creates an isolation layer of 
air above the free surface of the ferrogel. We 
have chosen the small aspect ratio F = /i/d ~ 5 
of the container, because the amount of ferrogel 
was limited. 

The vessel is flxed in the centre of a water 




(b) 

Figure 3: Experimental setup: (a) Scheme of 
the setup (not to scale); (b) photo of the tem- 
perature controlled container, flUed with ferro- 
gel. 
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cooled Helmholtz-pair-of coils (Oswald Mag- 
netteclinik Co.). An X-ray tube is mounted 
above the centre of the vessel at a distance of 
1606 mm. We measure the attenuation of X- 
rays passing through the ferrogel layer in the 
vertical direction with a photodiode array de- 
tector, which has a lateral resolution of 0.4 mm. 
Rays which are passing through crests are more 
attenuated than those passing through valleys 
(cf. Fig [3]a). The full surface topography is 
then reconstructed from the radioscopic images. 
The dynamic range of the detector of 16 bits 
translates into a vertical resolution of 0.4 |am. 
This resolution is achieved only when averaging 
many frames (~ 1000), because of a noisefloor 
with an RMS value of /irms = 20 \xm for a single 
frame. The calibration of the absolute height is 
limited to ~ 0.1mm by the stability of the X- 
ray source and the mechanical positioning. For 
details see Refs. [30,35]. 

Prior to each series of measurements, a batch 
of 12 ml of ferrogel is positioned in the empty 
vessel. This amount is molten by heating up the 
vessel to 90 °C, in this way creating a flat layer 
of ferrogel. Then the temperature is lowered to 
the desired value for the measurement and held 
constant. 



3 Results 

3.1 Quasistatic Experiments 

Figure H] presents characteristic topographies of 
the ferrogel for subcritical (a) and supercritical 
inductions (b). Due to the step-hke jump of 
the magnetisation at the container edge, a field 
gradient arises which attracts the ferrogel to- 
wards the container edge. Therefore a meniscus 
is formed. Chart (b) shows an example of the 
Rosensweig pattern. In the following, the height 
h of the central peak of the pattern serves as an 
order parameter and is estimated by fitting a 
paraboloid. 

Figure [5] gives the variation of the order param- 
eter for a slow increase (upward triangles) or 
decrease (downward triangles) of the magnetic 
induction for five different temperatures. In the 
interval from B = 19 mT to around 24 mT one 




Figure 4: A Rosensweig pattern emerging in fer- 
rogel, measured by radioscopy: (a) cross section 
for B = 22.9 mT; the transparent material illus- 
trates the size of the aluminium container; (b) 
full profile for B = 28.1 mT. The black con- 
tour lines are indicating consecutive levels with 
a distance of 1 mm. 
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notices a monotonic decay. In this regime no 
spikes exist. However, due to the formation of 
a meniscus (cf. Fig. Hla), the level of the ma- 
terial in the central part of the vessel, where h 
is estimated, decreases. For higher inductions, 
we observe a steep increase of h for all curves. 
It is this regime where the ferrogel spikes are 
emerging. 

Whereas in laterally infinitely extended liq- 
uid layers the transcritical bifurcation gives a 
proper scaling of the order parameter [22], for 
small pools of ferrofluid, imperfections induced 
by the container edges obscure the analytical 
scaling law. In this case, only a numerical model 
is available [36, 37]. In lack of an analytical 
expression, we can not extrapolate the values 
for the critical induction from a fit of h{B). 
As an approximation for the threshold of the 
predicted discontinuous transition we determine 
the induction -Bmax.up where the amplitude h{B) 
has its steepest inclination, i.e. dh/dB = max. 
It is determined from a spline fit and listed in 
table [T] for all investigated temperatures. A 
lower bound for the threshold is given analo- 
gously from the data for decreasing magnetic 
induction. It is denoted by -Bmax.dn- The hys- 
teresis, defined by ^^ax.up - ^max.dn, is in the 
range of 2 mT at 30 °C and shrinks to a frac- 
tion of 1 mT at 38 °C, as shown in the zoom 
presented in Fig. [5](b). 

Next we check, whether the time for a mea- 
surement cycle in B has an influence on the 
evolution of the order parameter. In Fig. [6] 
we present the results for a measurement cy- 
cle of 2 h (red dashed line) and of 4 h (blue 
solid line) for a temperature of 6^ = 30 °C, 
which is the lowest temperature, at which the 
measurements discussed beforehand have been 
performed. One clearly sees that the hystere- 
sis between the upward- and downward branch 
shrinks for longer cycle times. This variation is 
in the range of 1 mT. 

For higher temperatures (and lower G) the in- 
fluence of the cycle time becomes even less sig- 
nificant. Notwithstanding, the influence of the 
cycle time indicates, that we have different time 
scales in the ferrogel. These time scales will be 
studied next. 
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Figure 5: Height of the central peak for 
an increase (upward triangles) and decrease 
(downward triangles) of the magnetic induc- 
tion. The colour encodes the temperature from 
30 °C (blue) to 38 °C (red), as in the legend, 
(a) Full range. For clarity only every 10th data 
point is shown. The lines are splines to the full 
data set. The time for the measurement was 
2 h. (b) Zoom of the hysteresis. All data points 
are shown. 
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Figure 6: Height of the central peak for in- 
creasing induction (upward triangles) and for 
decreasing induction (downward triangles). The 
time between subsequent data points was dou- 
bled from 14.4 s (red dashed line) to 28.8 s (blue 
solid line). The time for a full cycle from OmT 
to 37 mT amounts to 2h and 4h, respectively. 



3.2 Dynamic Experiments 

3.2.1 Stress Relaxation Experiment 

In order to characterise the shear modulus G, 
we perform a stress relaxation experiment for a 
series of temperatures. We shear the sample by 
a deformation of 7 = 0.01 and measure the re- 
laxation of the stress r while the deformation is 
held constant. Figure [7] displays the data points 
recorded for various temperatures of the sample. 
For a specific temperature, the restoring force 
decays drastically during one second, which 
means that in our experiment, the shear modu- 
lus G = r/7 cannot be treated as constant. 
Commonly, a stretched exponential function is 
used to describe the time dependent moduli of 
linear viscoelastic media [38-42] 

Git,e) = Goi9)exp(^-(^-^yy (3) 

Here, the exponent (3 is restricted to the range 
[0, 1], with /3 = 1 for a simple exponential decay, 
and an increasingly broader distribution of re- 
laxation times for smaller values of /3. Moreover, 
this scaling law was recently observed in the 
stress relaxation of a triblock copolymer sub- 



Figure 7: Time dependent shear modulus 
G{t, 9) after a jump in the deformation of 7 = 
0.01 for different temperatures of the ferrogel. 
The temperatures have been varied from —10 °C 
to 32 °C. 



ject to an extensional strain [43], and consecu- 
tively explained in a model, which assumes that 
the copolymer reversibly splits into domains of 
different size [44]. These domains consist of 
a regular homogenous network of PS micelles 
which are interconnected by bridging chains of 
the middle block of the polymer. Also in our fer- 
rogels, we observe glassy PS micelles arranged 
in clusters, which are interconnected by bridg- 
ing chains of the ethylene-butylene middle block 
of the triblock copolymer gelator used. The size 
of these domains varies in the range from 60 to 
120 nm [6,33]. 

Next, we apply ([3]) to our relaxation data. 
The solid lines in Fig. [7] give approximations 
of G{t, 6) with (IH]) for different temperatures 
6 from 6 = 32 °C down to — 10°C. This tem- 
perature range was determined by the resolu- 
tion of our rheometer. All the fits use a com- 
mon exponent (3. The best value amounts to 
P = 0.34 ±0.01. 

The characteristic relaxation time to{6) of the 
ferrogel drops drastically with increasing tem- 
perature. It varies over six orders of magnitude, 
as shown in the Arrhenius plot Fig. [HI However, 
the dependence does not make up a straight line 
- therefore no simple Arrhenius behaviour with 
an activation energy can be inferred. 
Thus we observe for all temperatures a stretched 
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Figure 8: Arrhenius plot of the characteristic 
relaxation time. The data stem from fits of ([3]) 
to the measurements presented in Fig. [71 



exponential softening under load with a charac- 
teristic scaling exponent near Ys for all temper- 
atures, while the characteristic relaxation time 
to varies over six decades. These experimen- 
tal findings are in agreement with the obser- 
vation of domains [33] and the predicted scal- 
ing law based upon a reversible splitting of do- 
mains [44]. 

For t ^ oo, the time dependent shear modulus 
6) approaches 0. This means, that we have 
a viscoelastic soft magnetic material, without 
any long term elasticity. 

3.2.2 Periodic Driving Magnetically 

In the effort of constraining G to some finite 
value, we are next selecting a specific time scale. 
We periodically drive the imposed magnetic in- 
duction according to 



B{t) = Bo + AB sin(27r/z)t). 



(4) 



Here, Bq denotes the bias value, AB = 1.6 mT 
the driving amplitude, and /£> = 1 Hz the driv- 
ing frequency of the induction imposed by the 
Helmholtz-pair-of-coils. For each sample tem- 
perature, we measure the surface response for 
24 different values of Bq. For small Bq, the sur- 
face is oscillating with fu around its meniscus- 
like shape, which becomes alternatingly more 
and less pronounced concave. Beyond a thresh- 
old, spikes appear which are oscillating with the 



amplitude Ah around a mean height h with the 
driving frequency. For the whole range of val- 
ues, we observe a harmonic response which can 
be described by 



h{t) =h + Ahsm{2n fot + 



(5) 



In order to determine the quantities h and Ah, 
we measure the absorption of X-rays in the os- 
cillating surface pattern by means of an X-ray 
movie. For each data point we record 2200 
frames with a frame rate of 7.5 Hz. Each ab- 
sorption picture n is transformed into a height 
profile hn{x,y). From the series of height pro- 
files, we extract the desired quantities via the 
equations 



h{x,y) 

hcos{x,y) 
Ah{x,y) 

with 
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(6) 



The time delay between two consecutive frames 
amounts to At = (1/7.5) s, and W{t) = 
A^exp (— (t — ti/2)^/s^) denotes a normalised 
Gaussian window function with s = 0.4ti/2 and 
its centre at ti/2, i.e. at the half of the measured 
time interval. 

Figure M displays the time averaged height of 
the central extremum of the surface estimated 
in this way from the series of measurements at 
five different temperatures. In the interval from 
B = mT to around 24 mT one notices again 
a monotonic decrease of the central height, due 
to the growth of the meniscus at the container 
edge. For higher inductions, ferrogel spikes ap- 
pear which again lead to a steep increase of h 
for all curves. 

Also here, we use the maximal inclination of 
h{BQ), which is determined from a spline, as 
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Figure 9: Time averaged height h of the central 
peak of the surface pattern (cf. Fig. H]) versus 
the bias of the magnetic induction Bq. The hues 
stem from a sphne fit. The colour encodes the 
temperature as in Fig. O 



an estimate for the threshold. It is denoted 
by -Bmax.up and is listed in table [1] for all tem- 
peratures investigated. This estimate for the 
threshold is shifted for higher temperatures to- 
wards higher fields. The shift amounts only to 
^ 0.5 mT. 

Next, we present the oscillation amplitude 
of the central spike, as shown in Fig.[TO](a). The 
amplitude is increasing with increasing temper- 
atures, i.e. for softer ferrogel. From the max- 
imal increase, we determine an additional esti- 
mate for the threshold Ai?max.up- The results 
are listed in table [H Once more the shift for 
different temperatures is only ^ 0.5 mT. 
The X-ray movies of our oscillating spikes show, 
that the response is always harmonic and never 
sub-harmonic. This is also true for a driven 
harmonic oscillator. For strong damping and 
a driving frequency not much smaller than the 
resonant frequency, the amplitude of an oscilla- 
tor is inverse proportional to the damping con- 
stant. To uncover such a scaling in our measure- 
ments, we plot the rescaled amplitude Ah ■ \ri*\ 
in Fig.[Tn](b), where |?7*| = \G\/lj is the absolute 
value of the complex viscosity. Within the ex- 
perimental resolution, the graphs collapse onto 
a master curve. This indicates that the increase 
of Ah{B) under variation of 6 stems solely from 
the softening of the ferrogel. Therefore, one 
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(b) 

Figure 10: Oscillation amplitude versus the bias 
of the magnetic induction Bq. (a) Ah of the 
central peak (cf.Fig. Hj) (b) Rescaled amplitude 
Ah ■ \7]*\. The lines stem from a spline fit. The 
colour code for the temperature is the same as 
in Fig. [5] 
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possible simplified explanation for our experi- 
ment could be an oscillator, where the driving 
force comes from the magnetic stress at the edge 
of the container and the viscosity provides the 
damping. 

3.3 Comparison of the thresholds 
with predictions 

Now we want to take a closer look at the 
experimentally determined thresholds of the 
Rosensweig instability and how they compare 
to the predictions in Ref. [20]. For a ferrofluid, 
the threshold can be computed from the nonlin- 
ear magnetisation curve, the density p and the 
surface tension a via the linear stability analysis 
according to the book by Rosensweig [3], § 7.1. 
The critical magnetisation Mc of the fluid layer 
is given by 

Here = y/PchPt/po denotes the geometrical 
mean of the chord permeabihty = 
and the tangent permeability f^t = at the 

critical field. Together with M{H) and the 
jump condition of the magnetic field at the base 
of the dish, B = fiQ[H + M{H)], the critical in- 
duction can be determined from these implicit 
equations. 

For a ferrogel with a Hookean shear modulus 
G and a linear magnetisation curve, Ref. [20] 
provides the expression ([1]). We combine this 
equation with Eq. ([7]) for a non-linear M[H) to 
get the more general form 

Ml = -(l + -\{^+G). (8) 

If either one of the elasticity or non-linearity 
is left out, this equation reduces to ([1]) or ([7j), 
respectively. 

In table [H, we present the calculated critical 
inductions -BcFG according to Eq. (IHl). For 
all calculations, we utilise the storage modu- 
lus G"(lHz) as determined by the oscillatory 
measurements (cf. Fig. [1]). The other prop- 
erties entering Eq. ([8]) are taken from § 12.21 
accordingly. Specifically the surface tension is 



not well known for our gel. The value used 
here was ctff = 28.7 mN/m. Calculations for 
a = 35 mN/m show that a variation of a does 
not change -BcFG by more than 5 %. The mag- 
netisation curve has been modeled by the equa- 
tions given in [34]. 

The comparison between the experimentally de- 
termined thresholds in table [1] and the predic- 
tions in Ref. [20] reveals two prominent differ- 
ences. 

Firstly, we observe a decrease of the hysteresis 
of the pattern amplitude for a decrease of G, 
i.e. higher temperatures, when the magnetic in- 
duction is varied in a quasi-static manner. This 
hysteresis is denoted in Fig. [11] by the shaded 
and the upward and downward oriented full 
triangles. In contrast, Ref. [20] predicts an in- 
crease of the hysteresis under reduction of G for 
a ferrogel with Hookean elasticity. 
This can be explained from the temperature de- 
pendent relaxation process of the ferrogel. The 
relaxation times are increasing from r ^ 0.01s 
(at 32 °C) to T ^ 10000 s (at -10 °C). At the 
same time the measurement protocol for a full 
ramping of the magnetic induction was kept at 
2 h. The hysteresis increases because the mate- 
rial needs more and more time to follow a vari- 
ation of B. We have checked that by increasing 
the cycle time from 2h to 4h. The hysteresis 
was diminished, as shown in Fig. O 
Secondly, whereas G'(lHz) varies in the inves- 
tigated temperature range over two decades (cf. 
Fig.[T]), the threshold of the instability is vary- 
ing only within 10%. This is in contrast to the 
model [20] for a ferrogel with Hookean elastic- 
ity. However we found a different mechanical 
behaviour of our material. The shear modulus is 
strongly time- dependent and even vanishes for 

t — s> OO. 

In an attempt to select a specific, finite value for 
G, we modulated the magnetic induction with 
a fixed amplitude and a frequency of 1 Hz and 
applied this magnetic driving together with a 
magnetic bias induction Bq. Under increase of 
Bq a threshold is overcome and we observe a 
steep increase of the time averaged peak height 
h. The magnetic induction at which this transi- 
tion occurs is marked in Fig.[TT]by the blue filled 
circles. Also the oscillation amplitude Ah shows 



10 



25.5 F 




30 32 34 36 38 
Temperature 6 (°C) 



Figure 11: Temperature dependence of the 
thresholds. The black upward (downward) tri- 
angles mark -Bmax.up (-Bmax.dn) for an adiabatic 
increase (decrease), respectively. The shaded 
area in between denotes the hysteresis. The 
blue circles and red diamonds show -Bmax.up 
(ABjnax.up), respectively. The dashed line is just 
a guide for the eye. The theoretical value for a 
ferrofluid is displayed by the solid line. 




10 I i L.. i i i 1 

30 32 34 36 38 
Temperature 6 (°C) 

Figure 12: Comparison of theoretically (solid 
black line) and experimentally (filled circles) de- 
termined critical inductions, as listed in table [H 
The solid blue line is the theory for a ferrofluid. 
Note, that in this scale, all experimental values 
collapse onto that line. The dashed line denotes 
the vertical asymptote of the theoretical curve 
ate = 31.7 °C. 



a steep increase under variation of Bq. The cor- 
responding thresholds are shown in Fig. [11] by 
red diamonds. Within the experimental scatter 
both thresholds coincide. Similar to the static 
experiment, their values vary only in fractions of 
1 mT under decrease of the sample temperature 

e. 

In contrast, the theoretical value, estimated ac- 
cording to ([8]) at (j'(lHz), increases drastically 
under decrease of 6', as shown in FigJT^ Only 
for high temperatures, the experimental and 
theoretical values are close to each other. For 
lower temperatures, the gap between them in- 
creases drastically. According to ([8]) it even di- 
verges at a critical shear modulus of 

a = - = 77.8 Pa. (9) 

For our gel, the saturation magnetisation Ms 
of the ferrogel amounts to Ms = 14.7kA/m and 
the divergence occurs at ^ = 31.7 °C. Obviously, 
there is a mismatch of our complex, soft mate- 
rial and the linear model ([8]) when implementing 
G'(lHz). 

Let us look more closely on the small varia- 
tion of the experimentally determined thresh- 
olds in Fig. [m We see that under an increase 
of the temperature all estimates for the onset of 
the instability are slightly shifted to higher in- 
ductions. This shift can be understood from 
the fact that for increasing 6 the magnetisa- 
tion diminishes. Utilising the M{H, 6')-model 
by Ivanov et al. [34] , we have taken into account 
this effect. The corresponding results for -BcFF 
are marked in Fig. [11] by the solid line. The 
good agreement between the experiment and 
the plain Rosensweig estimate for the thresh- 
old computed in this way indicates that the 
shear modulus G'{1 Hz) has no influence on the 
threshold at all. 

4 Conclusion 

To conclude, we conducted the first measure- 
ments of the Rosensweig instability in thermore- 
versible ferrogels. The present material shows a 
complex viscoelastic relaxation process with an 
interesting critical exponent (3 ~ Ys, for exper- 
iments carried out in the time domain. It is 
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possibly explained by a reversible splitting of 
the polymer network into domains of different 
size [6,44]. Due to this relaxation, the threshold 
of the Rosensweig peaks is not much different 
from the Rosensweig instability in ferrofiuids 
under an adiabatic increase of a static magnetic 
field. However, the time scales are much slower. 
This is especially pronounced for the lower part 
of the investigated temperature range. Experi- 
ments with a periodic modulation of the mag- 
netic field show that the complex viscosity can 
be used to describe the response of the ferrogel. 
Such a complex elastic behaviour, however, is 
unsuitable to proof or rebut the model of the 
Rosensweig instability [20], which was derived 
for just a Hookean shear modulus. Certainly 
soft matter with finite G for t ^ oo would 
come closer to that model. One may think that 
a ferrogel with a higher gelator concentration 
will serve this aim. However, as a consequence 
for such a gel, the instability is completely hin- 
dered by the strong elastic modulus. Unfortu- 
nately, the amount of magnetite cannot be in- 
creased further to overcome the critical mag- 
netisation ([7]). Here, magnetic gels incorporat- 
ing particles with a higher saturation magneti- 
sation could solve this problem. Despite some 
attempts, with e.g. cobalt particles, such gels 
are not yet available. 
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Table 1: Critical inductions for various temperatures: -BcFF gives the results of for G = 
Pa, i.e. for a virtual ferrofluid with otherwise the material parameters of the ferrogel. -BcFG 
(mT) is estimated from G' at 1 Hz and the measured value of a. Experimentally, the thresholds 
were determined from the maximal incliniation of h{B), as marked by -Bmax.up for increasing B 
and -Bmax.dn for decreasing B. From the oscillatory measurements we determine the maximal 
inclination of h{B), marked by -Bmax.up? and the maximal inclination of Ah{B) abbreviated by 
A R 

'— ^-Djjiax.up 



e 


-BcFF 


-BcFG 


-Bmax.up 


-Bmax.dn 


-Bmax.up 


A -Bmax.up 


(°C) 


(mT) 


(mT) 


(mT) 


(mT) 


(mT) 


(mT) 


30 


24.76 


oo 


25.36 


23.43 


24.35 


24.40 


32 


24.90 


587.09 


24.99 


23.73 


24.36 


24.28 


34 


25.04 


83.17 


24.76 


24.02 


24.44 


24.39 


36 


25.18 


32.51 


24.81 


24.39 


24.64 


24.63 


38 


25.32 


28.03 


24.91 


24.62 


24.83 


24.87 
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1 Storage modulus G' (red circles) and loss modulus G" (blue squares) versus the 
temperature for the thermoreversible ferrogel 

2 Magnetisation as a function of the applied magnetic induction for the thermore- 
versible ferrogel. The solid line represents a least squares approximation with the 
model put forward by Ivanov [34]. The best fit parameters for the gamma distri- 
bution are: exponent a = 8.1 ± 1.3, mean diameter xq ■ {a + 1) = 8.1 ± 0.4 nm, 
volume fraction = 3.9 ± 0.09 %. The core magnetisation has been held constant 
at Mp = 446kA/m. The dashed line indicates the prediction for = 38 °C with 
the above parameters 

3 Experimental setup: (a) Scheme of the setup (not to scale); (b) photo of the 
temperature controlled container, filled with ferrogel 

4 A Rosensweig pattern emerging in ferrogel, measured by radioscopy: (a) cross sec- 
tion for B = 22.9 mT; the transparent material illustrates the size of the aluminium 
container; (b) full profile for B = 28.1 mT. The black contour lines are indicating 
consecutive levels with a distance of 1 mm 

5 Height of the central peak for an increase (upward triangles) and decrease (down- 
ward triangles) of the magnetic induction. The colour encodes the temperature 
from 30 °C (blue) to 38 °C (red), as in the legend, (a) Full range. For clarity only 
every 10th data point is shown. The lines are splines to the full data set. The time 
for the measurement was 2 h. (b) Zoom of the hysteresis. All data points are shown 

6 Height of the central peak for increasing induction (upward triangles) and for 
decreasing induction (downward triangles). The time between subsequent data 
points was doubled from 14.4s (red dashed line) to 28.8s (blue solid line). The 
time for a full cycle from OmT to 37 mT amounts to 2h and 4h, respectively. . . . 

7 Time dependent shear modulus G(t, 6) after a jump in the deformation of 7 = 0.01 
for different temperatures of the ferrogel. The temperatures have been varied from 
-10°C to 32°C 

8 Arrhenius plot of the characteristic relaxation time. The data stem from fits of ([2]) 
to the measurements presented in Fig. [71 

9 Time averaged height h of the central peak of the surface pattern (cf. Fig. H]) versus 
the bias of the magnetic induction Bq. The lines stem from a spline fit. The colour 
encodes the temperature as in Fig. [5] 

10 Oscillation amplitude versus the bias of the magnetic induction Bq. (a) A/i of the 
central peak (cf.Fig. H]) (b) Rescaled amplitude A/i ■ |?7*|. The lines stem from a 
spline fit. The colour code for the temperature is the same as in Fig. [S] 

11 Temperature dependence of the thresholds. The black upward (downward) trian- 
gles mark -Bmax.up (-Bmax.dn) for an adiabatic increase (decrease), respectively. The 
shaded area in between denotes the hysteresis. The blue circles and red diamonds 
show -Bmax.up (A-Bmax.up), respectively. The dashed line is just a guide for the eye. 
The theoretical value for a ferrofluid is displayed by the solid line 

12 Comparison of theoretically (solid black line) and experimentally (filled circles) 
determined critical inductions, as listed in tabled! The solid blue line is the theory 
for a ferrofluid. Note, that in this scale, all experimental values collapse onto that 
line. The dashed line denotes the vertical asymptote of the theoretical curve at 
e = 31.7°C 



16 



